Kinetic energy of uniform Bose-Einstein condensate 
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The Bogoliubov theory of a uniform weakly-nonideal Bose-Einstein condensate leads to a divergent 
expression for the kinetic energy of atoms. However, the latter can be determined provided that the 
dependence of the scattering length on atomic mass is known. The explicit expressions are derived 
for the kinetic energy through parameters that specify interatomic interaction. The kinetic energy 
of condensate atoms noticeably exceeds the total energy. Numerical data are presented for *^Rb 
and ^^Na condensates. 



The Bogoliubov theory of a weakly-nonideal (dilute) 
Bose-Einstein condensed gas is known to provide a quite 
good description of its properties |1| , permitting, in par- 
ticular, calculation of the total energy E of the gas in 
the ground state. However, the Bogoliubov approach 
can not be used for the calculation of the total kinetic 
energy, the important property of a condensate. Indeed, 
the distribution function of atoms out of condensate at 
large momenta p reads Np = (Airh^an) /p* where n is 
the gas density, and a is the scattering length of slow 
atoms (see, e.g. Sect. 4.3). 

Therefore, the total kinetic energy K , defined as the in- 



tegral K 
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diverges as J dp . This means 



that K is governed by small interatomic distances, of the 
order of the interatomic potential range tq, where the Bo- 
goliubov theory is inadequate. The total kinetic energy is 
an important quantity that characterizes the interaction, 
and in principle, it can be measured (see comments at 
the end of this letter). The calculation of this quantity is 
of considerable theoretical interest along with the total 
energy (note, that for an ideal Bose gas at T = 0, both 
E and K vanish). We emphasize however that we are 
speaking here about an uniform gas in which the kinetic 
energy is due to the presence of atoms outside the con- 
densate. This is a non-mean-field effect. An condensate 
in a trap possesses also kinetic energy due to its inhomo- 
geneity, which can be described by the mean-field theory. 

In order to calculate the kinetic energy integral, one 
should use a microscopic approach that explicitly takes 
into account properties of interatomic interaction. This 
can be done, for T = 0, by resorting to well-known for- 
mula that relates K to the derivative of E with respect 
to the mass of the atom, m, viz: K = —m (dE/dm) (see 
0, Problem to § 15). For the energy density (the energy 
per unit volume) one can use the first-order Bogoliubov 
expression E = 2ttTi'ii? (a/m) which is valid under con- 
dition na'^ ^ 1 (n is the density of atoms). In this way 
we get: 



n 2 d{a/m) 



K 

'e 



9 In {a/m) 



(1) 



dm ' E dm 
It is thus seen that K can be calculated provided the de- 



pendence of a on m is known. Notice that positiveness 
of the kinetic energy implies inequality on this depen- 
dence, d{a/m) I dm < 0. A simple example corresponds 
to atoms simulated by rigid spheres, when a does not de- 
pend on TO. In this case K = E, which is expected for 
this example. 

We begin from another simple, but instructive model. 
Let potential of interaction of atoms can be presented as 
a spherical square potential well of radius R and depth 
Uq . Then we have for the scattering length 



a = R- 



tan kR 



(2) 



where k = y/mU[)/m (see ^j, Problem 1 to § 132). Simple 
calculations give 



d {a/m) 



1 



3a + Rk^ {a - Rf 



dm 

For a shallow well, when kR <C 1, we find 
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and K 
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In this case <C ii^ • The non-zero value of the kinetic 
energy arises only in the second Born approximation. An 
interesting situation takes place when kR « (iV — 1/2) tt, 
where TV is a positive integer, and when the last term in 
the r.h.s. of Eq.(2) makes the main contribution to the 
scattering length. This is a case, when the well has a 
level with small binding energy, TV is the level number. 
Then a 3> i? and one gets from (P)): 
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For a more realistic case when the interaction between 
atoms is described by a smooth potential U (r) with a 
well, the mass dependence of a can be calculated within 
the quasiclassical approximation. The later is valid, when 
the number of levels N supported by the well is large 
and when the standard quasiclassical approximation for 
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bound states close to the dissociation threshold is cor- 
rected for a "nonstandard" matching of the WKB solu- 
tion in the well with a quantum solution in the region of 
the outer turning point Q. For our case, the potential 
in the region of outer turning is just a van der Waals 
potential ~a/r^ where a is the van der Waals constant. 
A relative corrections to this approximation are of the 
order of and we consider as a small parameter 
of our problem. For alkali dimers, N is of the order of 
100. 

The quasiclassical expression for the scattering length 
reads p: 

a = ao(l-tan$) (6) 
where ao is the "mean" scattering length 

a„ = 0.48(^)''' (7) 
and $ is the corrected quasiclassical phase integral 



$ - i / ^-2mU {r)dr - J . (8) 



Here is the inner turning point that corresponds to 
the classical motion with zero energy (the energy here is 
referred to the dissociation threshold). Under the quasi- 
classical conditions, the phase <I> is large and is related to 
N as ^ — ttN + {1). Eq.lO shows a very strong varia- 
tion of the scattering length upon a change of the phase 
integral by the increment of the order of unity. This is 
expected since the scattering length is a sensitive func- 
tion of a position of the upper, real or virtual, energy 
level supported by the potential well, see 0, § 133. In 
the calculation of the derivative da /dm we can therefore 
ignore a weak mass dependence of ao as given by Eq. jT)) 
and take into account the dependence of $ from Eq.© 
only. After simple derivations we get 
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Eq.® simplifies when the scattering length noticeably 
exceeds the "mean" scattering length oq, that is in the 
case of the resonance scattering on a weakly bound state. 
(Notice, that Bose-Einstein condensate is stable only if 
a > . Hence the level must be a real one.) For a ^ ag 
Eq.® assumes the form 



K = 



mao 



(10) 



The appearance of a large factor N in the expression 
for the kinetic energy, Eq. (|^- H10|I . reflects the accelera- 
tion of relative motion of a pair of atoms in the potential 
well. We note that the difference between iV^ and N be- 
haviour of K in Eqs. ^ and Hll|l comes from the fact 



that a square well and a realistic potentials show different 
pattern of convergence of the energy levels to the dissoci- 
ation limit. For any realistic potential, the energy levels 
converge to the limit in such a way that the spacing be- 
tween them decreases as \/—E, where E is the energy 
(negative) counted from the dissociation hmit 0- The 
square well does not show this behaviour, since even for 
high values of TV the bound state is not completely qua- 
siclassical because of non-quasiclassical behavior of the 
wave functions close to the walls. 

From Eqs.0 we get the following expression for the 
ratio of the kinetic to total energy: 



K _n Nao 
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This expression determines the K/E ratio via the scat- 
tering length a, "mean" scattering length ao and the 
number of bound states N. Two latter quantities, ao 
and N, characterize the microscopic parameters of the 
Bose-gas, in addition to the Bogoliubov "macroscopic" 
parameter a. Rough estimations of N for diatomic 
molecules can be done using the model of the Morse 
potential, N = 2D/huJe where D is dissociation energy 
and uje is the small-amplitude frequency of vibrations 
of a molecule. Using data from [3| we get TV = 135 
and N = 74: for *^Rb2 and ^^Na2 respectively. Ac- 
cepting these values for N, calculating ao from Q with 
a(Rb2)= 4770 a.u^^ a{Na2)= 1540 a.u. (see %, and 
borrowing a from [^j, U^, we get ao = 79as , a = lOOa^ 
for ®^Rb and ao = 43aB , a = 52a_B for ^"^Na. Finally we 
find K = 176E for ^^Rb and K = WOE for ^SNa. 

It follows from l|^- Hll|) that the kinetic energy in- 
creases with a faster than the total energy does, and K 
can be much more larger then E. (The mean-field ki- 
netic energy of inhomogeneous condensate is always less 
then the total energy.) In this case, the absolute value 
of the total potential energy, \V\ — \E — K\ , also ex- 
ceeds E , so that the total energy emerges from delicate 
balance of two large quantities which almost cancel each 
other. However, it should be taken into account that an 
increase in a, which happens when the last bound state 
approaches dissociation limit, is limited by the condition 
of a dilute gas, na'^ ^ 1. Under this condition, the total 
potential energy per atom is still much lower than the 
potential well depth. It is easy to show for the square 
well model. Indeed, for an dilute gas, one can conclude 
from © that /n - K/n - Uq (R/a) < L^o • Note 
that the kinetic energy actually always exceeds E, be- 
cause the potential energy is negative. Thus the positive 
total energy results at the expense of the kinetic energy. 

Large values of scattering amplitude have been reached 
in experiments near Feshbach resonances 0, 0| . How- 
ever in this case one has the "two-channel situation", 
when the resonance level is created by weak interaction 
between closed and open channels. This situation de- 
mands more careful theoretical considerations, which will 
be presented elsewhere. Notice only that the mass depen- 
dence of a for Rb isotopes was investigated numerically 
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in The kinetic energy of the condensate can be, 

in principle, measured. Indeed, according to the 
dynamical formfactor S {uj,q) of the condensate for a 
large enough momentum transfer q yields the momen- 
tum distribution; integration of this distribution allows 
one to calculate the kinetic energy. It might be more 
convenient to use the moments to„ (q) of the formfactor, 
™n (q) = / uj^ S {uj , q) duo . As shown in jl6l |. one can de- 
rive K from the asymptotic behavior (for g — s- cx)) of TO2 
or 7713. 

The dynamical formfactor of the condensate was first 
measured in by the two-photon Bragg spectroscopy. 
However, this method does not permit to get information 
on the momentum transfer of the order of Ti/a^ which is 
necessary to recover the kinetic energy. One can get in- 
formation on large values of the momentum transfer from 
neutron scattering data. In this way the kinetic energy 
of atoms in liquid *He was measured (sec, for example, 
[Tsf). However, this kind of experiments is difficult due to 
low density and small dimension of gas. More promising 
are experiments with scattering of fast atoms by conden- 
sates as discussed in Notwithstanding all difficulties, 



these types of experiments would open new important 
perspectives of investigation BEG at the atomic level. 

This theoretical approach can be also compared with 
results of the Quantum Monte Carlo calculations. First 
calculations of the kinetic energy were presented in [l^ . 
The authors used a model with the delta-potential inter- 
action and kinetic energy they calculated was related to 
the gas inhomogeneity in a trap. However, it is possible 
to perform such calculations for a realistic potential. 

In conclusion, we have calculated the kinetic energy of 
a uniform Bose-Einstein condensate in its ground state. 
This quantity, which diverges in the Bogoliubov theory, 
is determined by the microscopic properties of a dilute 
condensed gas: the dependence of the scattering length 
on the atomic mass. Different possibility of measurement 
the kinetic energy of the condensate are discussed. 
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